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Abstract. We confront our quasi-particle model for the equation of state of strongly interacting matter 
with recent first-principle QCD calculations. In particular, we test its applicability at finite baryon densities 
by comparing with Taylor expansion coefficients of the pressure for two quark flavours. We outline a chain 
of approximations starting from the ^-functional approach to QCD which motivates the quasi-particle 
picture. 



1 Introduction 

In the last years, great progress has been made in the 
numerical evaluation of QCD thermodynamics from first 
principles (dubbed lattice QCD) even for finite chemical 
potentials PEEEE • While various perturbative expan- 
sions 6,7,8,9,10 fail in describing thermodynamics of 
strongly interacting matter in the vicinity of T c (T c be- 



ing the (pseudo-) critical temperature of deconfinement 
, and chiral symmetry restauration), different phenomeno- 
logical approaches exist which aim to reproduce the non- 
\ perturbative behaviour. For instance, models based on 
quasi-particle pictures with effectively modified proper- 
ties due to strong interactions are successful in describing 
lattice QCD results 1111112111311141115111611171 . Analytical ap- 
proaches with a rigorous link to QCD (cf. [TBI for a survey) 
such as direct HTL resummation 19 20 or ^-functional 
approach [2Tll22U23j formulated in terms of dressed prop- 
agators are successful in describing lattice QCD on tem- 
peratures T > 2T C . 

It is the aim of the present paper to show the success- 
ful applicability of our quasi-particle model (QPM) for 
describing lattice QCD results and to motivate the model 
starting from the ^-functional approach to QCD. In sec- 
tion[2J we review the QPM and compare with recent lattice 
QCD results for pressure and entropy density. In section 
03 a possible chain of approximations is outlined starting 
from QCD within the ^-functional approximation scheme 
which motivates our formulation of QCD thermodynamics 
in terms of quasi-particle excitations. We summarize our 
results in section 2] 



2 QPM and Comparison with Lattice QCD 

In our model, the pressure p for Nf = 2 light quark 
flavours in thermal equilibrium as function of tempera- 



ture T and one chemical potential \i q (fig = 0) reads 

p(T, H q ) = ^2 Pa- B(T, H q ) , 



(1) 



where p a = d a /(6ir 2 ) / °° dkk 4 (/+ + f~) /u a denote the 
partial pressures of quarks (q) and transverse gluons (g). 
Here, d q = 2N f N c , d g = N 2 - 1, N c = 3, and /± = 
(exp([cj a =p l^a]/T) + 5a) -1 with S q = 1 for fermions and 
S g — —1 for bosons. B(T,/j, q ) is determined from ther- 
modynamic self-consistency and the stationarity of p un- 
der functional variation with respect to the self-energies, 
Sp/5II a = . II a enter the quasi-particle dispersion re- 
lations uo a being approximated by asymptotic mass shell 
expressions near the light cone, uo a = \/k 2 + II a . We em- 
ploy the asymptotic expressions of the gauge independent 
hard thermal (dense) loop self-energies [2SJ- Finite bare 
quark masses mo- q as used in lattice simulations can be 
implemented following 

By replacing the running coupling g 2 in U a with an ef- 
fective coupling G 2 (T, \i q \ non-perturbative effects in the 
vicinity of T c are accomodated. In this way, we achieve 
enough flexibility to describe lattice QCD results. We 
parametrize G 2 (T 1 pL q — 0) (cf. 27 for details) by 



G 2 (T, N =0) 



G 2 (£(T)), T>T C , 



r (2)V 
r (2)( 



G 2 (£(T C ))+6(1 



rp ), T <C Tq 



(2) 



where G 2 2 ^ is the relevant part of the 2-loop running cou- 
pling and £(T) = A(T — T s )/T c contains a scale param- 
eter A and an infrared regulator T s . The effective cou- 
pling G 2 for arbitrary T and fi q can be found by solving 
a quasi-linear partial differential equation which follows 
from Maxwell's relation, 



dG 2 

dfl q 



dG 2 1 



(3) 
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The coefficients in Q explicitly read (neglecting for sim- 
plicity additional contributions stemming from In- 
dependent bare quark masses as employed in lattice sim- 
ulations) 



(4) 
(5) 




N C N S 



(6) 



Here, 



(/,-)■' -^ IT ui? i,(7) 



L q ft 



e « 



IT 



J a 



1 e<l T 



h = 



2T 
d 



h 2 / T 



1 



(8) 



(9) 



with ff = f g , uj± 



± fi q , L a = 2k 2 /3 + n a /2 and 



Cf = (N* - 1)/2N C . 

Entropy density s = dp/dT = ^ a=q g s a and net den- 
sity n = n q = dp/dfi q follow from as 



4- TT 



2<k 2 T J v 

-Ma [fa ~ fa] ) > 



q ~ 2tT 2 



^ y o dkk 2 [/+-/-] . 



[/+ + /"] (10) 



(ii) 



In Fig. d we exhibit QPM results for p and s at /i = 
compared with lattice QCD results for different numbers 
of quark flavours 28,29). 

Recently, the decomposition of p into a Taylor series in 
powers of (ji q /T) for small \i q was studied in lattice QCD 

inn, 

CO 

p(T, M9 )=T 4 £c„(T)(^) . (12) 



77=0 



The expansion coefficients c„(T), vanishing for odd n and 
depending only on temperature T, follow using |JTJ from 



c n (T) = 



1 d n {p/T A 



n\ d(fi q /T) n 



(13) 





Fig. 1. Comparison of our QPM with lattice QCD results 
(symbols) for p(T,p q — 0)/T 4 (upper panel) and s(T,p q — 
0)/T 3 (lower panel) as functions of T/T c for Nf = 2 (squares) 
EH] and N f = 2 + 1 (p s = 0) (circles) Raw lattice QCD 

data are continuum extrapolated as advocated in 28 30 . QPM 
parameters: A = 4.4, T s = 0.67T C , b = 344.4, B(T C ) = 0.31T C 4 
with T c 175 MeV for N f = 2 and A = 7.6, T s 0.80T C , 
6 = 348.7, B(T C ) = 0.52T C 4 with T c = 170 MeV for N f = 2 + 1. 



c n (T) depend on G 2 and its derivatives with respect to \i q 
at /i g = 0, thus testing (J3J). Furthermore, net density n can 
also be decomposed into a Taylor series at small \i q with 
expansion coefficients c n (T). Therefore, the higher order 
coefficients C2,a$(T) serve for a more direct test of the ap- 
plicability of our model at finite \i q . In Fig. [21 we compare 
C2,4,6(^) evaluated from (p"3|) with lattice QCD results for 
Nf = 2 jSI]. In particular, the pronounced structures in 
the vicinity of T c are fairly well reproduced (cf. |27|). 



3 Foundations of the QPM 

Having successfully reproduced first-principle lattice QCD 
results, it would be desirable to establish contact between 
our ad hoc introduced QPM in section [21 and QCD as the 
fundamental microscopic gauge field theory of strong in- 
teractions. In order to motivate our quasi-particle model, 
we present a possible chain of approximations starting 
from QCD within the ^-functional approach following the 
pioneering work 21,22,23 . We concentrate on entropy 
density s and net density n, as they turn out to possess a 
simple structure supporting the picture of quasi-particle 
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Fig. 2. Comparison of our QPM with lattice QCD results 
(symbols) (31) for C2,a(T) (upper panel) and cq(T) (lower 
panel) as functions of T/T c for Nf = 2. QPM parameters: 
A = 12.0, T s = 0.87T C , b = 426.1, with T c = 175 MeV. The 
horizontal lines at T > T c depict the corresponding Stefan- 
Boltzmann values highlighting the effects of strong interaction 
near T c . 



excitations. Other thermodynamic quantities such as pres- 
sure p or energy density e are determined from s and n. 
Although rather strong assumptions become mandatory 
in the derivation, one should be aware of the remarkable 
success of our QPM in describing lattice QCD results. 

In the ^-functional approach pSIM| to QCD, the ther- 
modynamic potential Q — —T In Z can be expressed as a 
functional of dressed propagators of gluons D, quarks S 
and Faddeev-Popov ghost fields G, 



Q[D,S, G] 



^TrpnD- 1 - 
-TrpnG- 1 



nD]-Tr[\iiS- 1 -US] (14) 
-EG\+$[D,S, G] . 



Here, ghost field contributions compensate for possible un- 
physical degrees of freedom in the gluon propagator. While 
the propagators in (jTljl depend on the specific gauge, 
Q = —pV must be gauge independent. For convenience, 
we choose the Coulomb gauge in the following in which 
ghost fields do not propagate and the gluon propagator 
consists only of the physical transverse and longitudinal 
modes. The functional ^[D, S] is given by the infinite sum 
of all 2 -particle irreducible skeleton diagrams constructed 
from D and S. 



The self-energies are related to the dressed propagators 
by Dyson's equations 



n[D) = D' 1 - Dq 1 , U[S] = S' 1 - s 



(15) 



where Do and So represent the bare propagators of gluon 
and quark fields, respectively. Demanding the stationarity 
of Q under functional variation with respect to the dressed 
propagators j^U 



5f2[D, S] 



5D 



6fi[D,S\ 



Do 



ss 



o, 



(16) 



the self-energies follow self-consistently by cutting a dressed 
propagator line in ^ resulting in the gap equations 

n _ a «3M >E __«IM. (17) 

The trace "Tr" in ((14)) has to be taken over all states 
of the relativistic many-particle system. In the imaginary 
time formalism it can be rewritten in the form Tr — > 
trPVTY,n=-oc Jd 3 k/(27r) 3 . Here, V is the volume of the 
system, (3 = 1/T and "tr" denotes the remaining trace 
over occuring discrete indices including colour, flavour, 
Lorentz or spinor indices. Introducing the four-momentum 
k v = (u;,k) = {iuo ni k), the sums have to be taken over the 
Matsubara frequencies uo n = 2uttT (or (2n + 1)ttT — iji) 
for gluons (or quarks). They can be evaluated by using 
standard contour integration techniques in the complex 
cj-plane 25 ,35 wrapping up the poles of the propagators. 
Expressing the analytic propagators in terms of their spec- 
tral densities p, one can define 



Pd(s)(v, \k\) 



21im Im D(S)(uj 



■ie,\k\) (18) 



for real uj. Similarly, the imaginary parts of functions of 
the analytic propagators obeying the same pole structures 
can be defined. Hence, £2 reads with retarded propagators 
D and S depending on uo and k = |fe| 



tr/ 



d 4 fe 



+2tTf^f(u)Tm\hiS 



(^)Im[lnD- 1 -IID] (19) 



2S] + £$[D,S\. 



where / d A k = j d 3 k j dw, andn(cj) = (e^-l)" 1 (f(u) = 
( 6 /30-m) _|_ 1) _1 ) denotes the statistical distribution func- 
tion for gluons (quarks with chemical potential pi). 

Due to the stationarity property (|16|). entropy den- 
sity s = -d{Q/V)/dT and net density n = -d(f2/V)/dfi 
contain only explicit temperature and chemical potential 
derivatives of n{uj) and /(a;), although the propagators in 
(|T9j) depend implicitly on T and pi through their spectral 
densities. Using Im (IID) = Im II Re D + ReTIImD, one 
finds for the entropy density s = s g + s q + s' with 



-tr 



-2tr 



d A k dn(<jj) 
(2tt) 4 dT 



[ImlnD- 1 -ImllReD} , (20) 



d A k df(u) 
(2tt) 4 dT 



[ImlnS- 1 -Im^Re^] , (21) 
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g(ig[Ag]) 



dT 



D,S 



2tr/ 



d 4 k dfju) 
(2tt) 4 8T 



(22) 



Similarly, for the net density one finds n = n q + n' with 
d 4 k df(u) 



n q = — 2tr 
n 2tr 



(2tt) 4 <9/i 
d 4 /c dfjuj) 
(2tt) 4 <9/x 



[IinkiS -1 -ImZ , Re5'],(23) 
Re^ImS (24) 



D,5 



While the sum integrals in ft (p"9|) contain ultraviolet di- 
vergencies which must be regularized, the expressions for 
s g , s q and n q in 12T)) and ((2HJ) are manifestly ultravi- 
olet convergent because the derivatives of the statistical 
distribution functions vanish for uj — > ±oo. In addition, 
introducing real multiplicative renormalization factors for 
propagators and self-energies, these factors simply drop 
out of s and n. 

Self-consistent (or ^-derivable) approximation schemes 
preserve the stationarity property (|16[) of ft when truncat- 
ing the infinite sum in ^ at a specific loop order while cor- 
responding self-energies and propagators are self- 
consistently evaluated from ([T7)l and Dyson's equations. 
Nevertheless, self-consistency does not guarantee gauge 
invariance which is an important issue in truncated ex- 
pansion schemes. In fact, by modifying propagators but 
leaving vertices unaffected Ward identities are violated. 

We consider @ at 2 - loop order in the following which 
is diagrammatically represented by jSHJ 



<2> 



(25) 



Here, wiggly (solid) lines denote gluons (quarks). The self- 
consistent self-energies are accordingly 



n 




(26) 



(27) 



Although vertex corrections can be implemented self- 
consistently [37| , they turn out to be negligible at 2 - loop 
order in ^ [52 . In addition, s' = n' = is found for the 
residual contributions of entropy density and net density 
in {23 and (EH) at 2 - loop order 22 j. This topological 
feature being related to (JT7jl has also been observed in 
massless <£ 4 -theory and in QED [3HJ. 



Concentrating on the gluonic contribution s g , ((SDJ) can 
be rewritten by using the identity 

LnpnD- 1 ^,*;)] = -Trsgn(^) ©(-RelT 1 ^ k)) (28) 

/ Im77(c<;, k) 



- arctan 



yReD- 1 ^, k) 



where — tt/2 < arctan x < tt/2. Hence, s g can be decom- 
posed into s g = s Qy Qp + Sg^LD with 



S 9,QP 



s q ,ld 



duo dn(uj] 
~2 df~ 



r d 3 k r 
J (2^F7 

f d A k dn(uo) r 
tr / - — — ^ v J UmllReD 
J (2^) 4 dT \ 

( imn ^ I 



■sgn(cj)6)(-ReL>- 1 ),(29) 



(30) 



- arctan 



VReL>- 



Here, (|29|) accounts for the contribution of dynamical 
quasi-particles to s g defined by the poles of D and 
represents the contribution from the continuum part of 
the spectral density associated with a cut below the light 
cone | a; | < k 26 , 40 representing Landau damping. Ap- 
plying a similar identity for Im[ln S^ 1 ^, fc)], s q and n q in 
(|2ll l23|) can be decomposed similarly into quasi-particle 
and Landau-damping contributions. 

In Coulomb gauge, D consists of a longitudinal and 
a transverse part, Dl and Dt- Similarly, the (massless) 
quark propagator consists of two different branches with 
chirality either equal (positive energy states) or oppo- 
site (negative energy states) to helicity. By employing the 
gauge invariant hard thermal loop (HTL) expressions 77 
(XJ) for the gluon (quark) self-energies in the following, 
one obtains gauge invariant approximations of s and n. 
The HTL expressions read 25 





= m 2 D 




n T (co, k) 


-K 


m 2 D + 




_ M 2 
k 


0- 



u> , u> + k 



2k uj — k J ' 
to 2 - k 2 



(31) 



k 2 



n L (u,k)\, (32) 



uj =p k _ + /c 
— — — In - 

2/c uj — k 



(33) 



with Debye screening mass (allowing, in general, for dif- 
ferent chemical potentials jii) 



[2N c + N f ]T 2 + N c J2^] (34) 



long- wavelength fermionic frequency 
JV?-1 



M 



16iV, 



5 



(35) 



and running coupling g 2 . Although being derived origi- 
nally for soft external momenta uj, k ~ gT <C T, they 
coincide on the light cone with complete 1-loop results 
0T] as exhibited in Fig. 03 Finite quark masses, m < T, 
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modes, one finds 



Fig. 3. Comparison of HTL approximation (the nearly indis- 
tinguishable short and long dashed curves are for massless (pH)) 
(cf. gH) and massive quarks, m q — 0.4 T, as used in lattice cal- 
culations 28 29 30 31 , respectively) with 1-loop results (solid 
curves for massive quarks) |42| of the scaled real part of the 
transverse gluon self-energy as function of ou/T for k/T = 0.1, 
0.5 and 1 from left to right. Nf = 2, {ii = 0. In the vicinity 
of the light cone uu = k (shaded regions), HTL results ap- 
proximate 1-loop results fairly well even for k ~ T which is 
the thermodynamically relevant part in momentum space but 
with decreasing agreement for increasing k. 



turn out to be negligible. The corresponding propagators 
are evaluated from Dyson's equations. 

For k ~ T, ii the poles of both, longitudinal gluon 
propagator as well as abnormal fermion branch have ex- 
ponentially vanishing residues 26 giving only minor con- 
tributions to the thermodynamics. Therefore, we assume 
that these collective modes can be neglected in the fol- 
lowing. Furthermore, being a severe approximation, we 
also neglect any imaginary parts of the self-energies, i. e. 
IulZTt = ImZ+ = 0. Then, the Landau damping contribu- 
tions to s g , s q and n q vanish. Finite width effects associ- 
ated with imaginary parts of the self-energies are discussed 
by Peshier HU|. Including Landau damping as well as the 
exponentially suppressed modes, it was shown in [H] that 
in this way some ambiguities arising when solving (0 can 
be eliminated. 



Performing the c^-integration in ((29)) (but now for Dt), 
the only contributions stem from uu 2 > uu\ because of the 
(9-function, where uut is the positive solution of uu 2 — k 2 — 
77t(cj, k) = 0. Therefore, the cj-integral in (|2HJ) reads 



dujdn(uj) * ! 

"2 — 7^sgn(^)(9(-ReL> T ) 



dn(uu) 



(36) 



Joo "2"V 



<^ T dto fdn(-u) 
2 V dT 



dT 



, Qp = _ 2(JV * _ i)| (i n(1 _ e -^) (37 ) 



(3ujt 



of3u T _ 2 



Similarly, s Qj qp can be evaluated, where non- vanishing 
contributions to the cj-integration stem from uu > uu+. 
Here, uu+ is the solution of uu — k — £+(cc;, k) = for the 
positive fermion branch. Using —df(uu)/dT = da(uu)/duu 
for the spectral function a(uu) = — f(uu) In f(uu) — (1 — 
/(o;))ln(l — f(uu)), the a;-integral can be integrated by 
parts. Antiquarks are included by simply replacing \i — > 
—fi in f(uu). Taking the trace over remaining spin, colour 
and flavour degrees of freedom, one finds 



S q,QP 



2N c N f J A ( ln(l + e -/^+-M)) (38) 



(27T) 

(3(u) + -fj,) 



+2N c N f J ^ ^ In(l + e-^+^ ) 



/3(c^++/i) 



s g,QP an d Sq,QP in (|37ll55j) represent entropy density con- 
tributions of non-interacting quasi-particles with quantum 
numbers of transverse gluons (quarks) and dispersion re- 
lation UUt (^+)- 

Correspondingly, nq,QP is evaluated using 
—df(uu)/dji = df(uu)/duu. Adding antiquarks by \i — > — /i 
in f(uu) (note, now df(uu)/dfi = df(uu)/duu) and taking 
the trace, n q ^Qp reads 



n qiQP = 2N c N f J 



d 6 k 
(2tt)3 



+ 1 



(39) 



Finally, we approximate the quasi-particle dispersion re- 
lations by the asymptotic mass shell expressions near the 
light cone, thus neglecting any momentum or energy de- 
pendence of the self-energies. We employ uut — > uo g = 
xjk 2 + n g and uu+ — > uu q = \Jk 2 + II q as in section El with 
asymptotic masses II g = rh 2 D /2 and II q = 2M 2 , consider- 
ing only one chemical potential = /i g (^ /if — > Nffi 2 ). 
Integrating the logarithmic terms in (|37l 135)) by parts 
(note that ((39)) already obeys the desired form), one ex- 
actly recovers expressions ([TU1ITT]) of our QPM, where the 
replacement of g 2 by G 2 remains as phenomenological pro- 
cedure on top of the listed " approximations" . 



The remaining integration is performed through an in- 
tegration by parts using —dn(uu)/dT = dn(—uu)/dT = 
dcr(uu) I duu for the spectral function <j(uu) = —n(uu) In n(uu) + 
(1 + n(uu)) In (1 + n{uu)). Taking the trace over polariza- 
tion and colour degrees of freedom for the transverse gluon 



4 Conclusions 

In summary, motivated by the successful reproduction of 
available results of QCD thermodynamics, we attempted 
a collection of necessary steps to establish the link of our 
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employed quasi-particle model to QCD. Quite severe as- 
sumptions had to be made. Even with these, resulting in 
the formal structure of our model, an additional and cru- 
cial point is the parametrization of the effective coupling. 
While allowing for an accurate two-parameter fit of many 
different lattice QCD results, it requires a foundation. In 
this respect, we refer to the work in 43 , where the authors 
argue that the pure quasi-particle excitations, deduced 
from a preliminary study of the poles of quark and gluon 
propagators 03j, are too heavy to saturate the pressure 
delivered from lattice calculations, e. g. 29 , i. e. signalling 
the necessity of including additional degrees of freedom. 
Further systematic studies of the relevant degrees of free- 
dom in the strongly coupled quark-gluon medium near T c 
are highly awaited to have some guidance. 

An additional issue is the chiral extrapolation. The 
quark masses of 0.4 T employed in the here analyzed lat- 
tice simulations correspond to unphysically heavy pions 
with ~ 770 MeV. In the 1-loop and HTL gluon self- 
energy considered here, such a finite quark mass has a 
tiny (negligible) impact. A more rigorous treatment of fi- 
nite quark mass effects must be accomplished to arrive at 
a suitable chiral extrapolation. 
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